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Abstract 

We report measurements of the Coulomb-blockade current in resistively 
isolated {Risoi ^ h/e^) tunnel junctions for the temperature range 60 mK 
< T < 230 mK where the charging energy is much greater than the ther- 
mal energy. A zero-bias resistance i?o of up to 10^i?T (the tunnel resistance 
of the bare junction) is obtained. For eV <C Ec, the I-V curves for a given 
Risoi scale as a function of V/T, with / oc y^(^^soi) over a range of V. The 
data agree well with numerical calculations of the tunneling rate that include 
environmental effects. 

PACS number: 73.23.Hk, 73,40,Rw, 71.10.Pm, 73.40. Gk 
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Coulomb blockade in small-capacitance tunnel junctions has been intensively studied in 
the past decade In general, Coulomb blockade is observed when the electric charge Q 

is localized on a small conducting island, i.e. when both thermal and quantum fluctuations 
of the charge are suppressed. Thermal fluctuations are suppressed when the temperature is 
below the charging energy Ec = e^/2C of the island capacitance C. Suppression of quantum 
fluctuations requires the isolation resistance Risoi between the island and the "outside world" 
to be much greater than the quantum resistance Rk = h/e^. In systems of several tunnel 
junctions, where the electrodes are physically isolated by tunnel barriers, this condition 
is easily satisfled, leading to strong Coulomb blockade 0]. However, in a single, current- 
biased tunnel junction, the localization of charge can only be achieved in practice with an 
Ohmic resistor [Q. This requires the fabrication of a compact resistor of large magnitude, 
Risoi ^ Rk- The major problem in achieving this is that the isolating resistor acts as a 
lossy (RC) transmission line (unless it is very short), with an impedance Zjsoiiyj) ^ Rk for 
frequencies approaching Ec/h, which are important for the tunneling process. Previously 
published reports on resistively isolated junctions p-pi?| showed relatively small increases in 
the zero-bias resistance Rq compared to the high-bias tunnel resistance Rt- In this Letter we 
report the observation of strong {Rq/Rt up to 10^) Coulomb blockade in a resistively isolated 
junction, and the flrst observations of the predicted power-law variation of the tunneling 
current I{V) in the blockade region with an exponent dependent on Rjsoh iii analogy to 
recent observations of similar behavior in Luttinger liquids [pJ]-[T^. This strong blockade 
enables us to quantitatively test the "environment" theory, based on the electrodynamic 
description of the resistor. 

The theory for the effect of such resistive environments (characterized by continuous 
charge transfer) on the tunneling rate through the junction is now well developed for the 
case Rt ^ Rk p6Hl9[] . These calculations predict that at zero temperature and for V <^Vc 
{Vc = Ec/e) the current varies as a power of energy change E of the system when an electron 
tunnels. Since, for a single junction, E = eV, this gives: 
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J oc oc y", a = + 1. (1) 

Rk 

For finite (but sufficiently low) temperatures, the zero-bias resistance Rq varies as a power 
T : 

Ro oc T^-". (2) 

The origin of these power laws is the one- dimensional nature, i.e. the linear variation 
with frequency, of the number of the "photon modes" used to model the resistors. In this 
respect, they are analogous to the power law (/ oc V°^^) recently predicted and observed for 
tunneling in Luttinger liquids ||Tl|-[T4|. The two systems are very different on the microscopic 



level: a strictly ID ballistic conductor for the Luttinger liquid, and a disordered macroscopic 
conductor with a large number of transverse electron modes leading to 3D electron motion 
for the Ohmic resistor. Nonetheless, the correlation functions of both systems are dominated 
by bosonic excitations with a linear spectrum. This particular spectrum gives rise to the 
power-law I-V characteristics. 

The sample (see upper inset of Fig.|l]) is fabricated with a three-level process. First, 
the resistor level is patterned and a bilayer film of 3nm Cr and 25nm of Au is sequentially 
deposited without breaking vacuum. An etch window is then patterned over parts of the 
leads and Au removed to form the Cr wires, which are approximately lOOnm wide and 
have a resistivity of about 10 kQ/fim. Finally, the junction- level mask is patterned, and the 
standard shadow-evaporation technique with a vertical offset is used to form the Al/AlO^^/Al 
tunnel junctions (dark rectangles in upper inset). The specific capacitance of 4.5 fF //im^ 
gives a nominal capacitance for the small junction of about C = 230 aF. Adding to this 
is the parasitic capacitance Cp of the leads between the junction and the resistors, which 
we estimate to be ~35 aF based on calculations of the self capacitance of other structures 
with similar wires The spaces between the resistors are 1.5 fim horizontally and 1.5 
fim vertically. Calculations show that the isolating resistors can be treated as ideal lumped 
resistors for frequencies up to 100 GHz. The resistors have nearly temperature- independent. 
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linear I-V characteristics to the lowest temperature. Their noise properties, measured down 
to 1.2 K, are found to be consistent with Johnson noise for the range of bias currents 
used in the data presented below. Measurements are made with the sample enclosed in 
a liquid-He-filled copper cell located on the temperature-regulated platform of a dilution 
refrigerator. All leads entering this cell pass though high-attenuation microwave filters ^1 



which are thermally anchored to the platform. A magnetic field of 1 T is applied to suppress 
superconductivity in the sample. 

The fabrication process, chosen to minimize Cp, results in a secondary series junction with 
an area (and capacitance Cs) about 15 times that of the small junction, making the device 
into an asymmetric transistor, albeit one whose conductance is almost entirely determined 
by the small junction. The energy change as an electron tunnels through the small junction is 
given by E{q) = hi[eV + {q — e/2) /Cs], where q is the island charge and k = Cg/ {C + Cg) ~ 1. 
For the region of interest, V <^Vc, where the current is determined by the forward tunneling 
rate through the small junction, the voltage for a fixed bias current is minimized for q = e/2. 
In this case, E{e/2) = neV and, within the sequential tunneling model, the asymmetric 
transistor is predicted to behave essentially as a single junction |TP|: It has no blockade for 



Risoi 0, has a blockade voltage = e/2C when Rjsoi ^ Rk, and exhibits the power-law 
behavior of Eqs.|l] and ^ with a K^a. Both the smaller "bare" resistance of the secondary 
junction and the much weaker effect of the high impedance environment (k^ < 0.01) ensure 
that the first junction is the bottleneck for the electron transfer across the device. If the 
resistance Rg of the (probably leaky) secondary junction is much less than h/e"^, k, ^ 1, 
and the correspondence between the transistor with q = e/2 and the single junction is even 
closer, with Rg simply increasing Risoi slightly. 

We present data from two samples with nominally identical junctions but with isolating 
resistors of different length and resistance: 4 /im (40kf2) for sample 1 and 8 fim (75 kfl) 
for sample 2. The parameters for these samples are summarized in Table |I[ All data are 
taken with q = e/2 as determined by adjusting the back-plane gate to minimize V for 
V Vc- A low-temperature I-V curve for sample 1 is shown in Fig.|l]. The voltage offset. 



Vc = 455 /iV, due to Coulomb blockade is clearly seen in the linear-scale lower inset. The 
current variation for \^ < is more clearly seen in the main log-log plot, where it is compared 
with the numerical solution of the theory. The large current suppression {Rq > 4 x 10'^ -Rt) 
indicates the strong blockade seen in this sample. Figure shows I-V curves taken at several 
temperatures for both samples, again compared with the result of numerical solution. One 
can see that sample 2, which has a larger Risoi, has a stronger blockade and a steeper 
power-law region at the lowest temperature. 

The current I(y, T) (dashed lines in Figs. and |]) is calculated using the standard golden- 
rule approach with the effects of the environment included through the correlation function 
of phase fluctuations due to the total environmental impedance Zt{uj) [p!^ -[T8|. We take the 
real part of this impedance to be: Re[Zt(c<j)] = Risoi/ (1 + i'-^/^B.c)'^) where ujic = ^/RisoiCe 
and Ce is the effective junction capacitance corresponding to the experimentally determined 
value of Ec (see Table 1). Ec and Rt are determined from the high- voltage data, where the 
corrections due to finite temperature and Risoi are negligible. The isolating resistance Risoi is 
then used as the only adjustable parameter in fitting the data for V < Vc, giving an excellent 
fit over six orders of magnitude in the current. In view of the essentially identical behavior 
of the transistor with q = e/2 and the single junction, and the impossibility of accurately 
determining Rg, we treat the sample as a single junction and admit some uncertainty in 
our comparison of the measured and fitted values of Risoi- Table H compares this best-fit 
value of Risoi with the independently measured value. For both samples the best-fit value is 
approximately 10 kQ higher than the measured value. This may very well be the resistance 
of R, 

As noted, the physical basis for the power-law behaviors seen in these samples and in 
Luttinger liquid is the same. Therefore it is of considerable interest to examine in more detail 
the region of validity for those behaviors in our devices. While the full theory (Figs. ^ and 0) 
requires numerical integration, analytical expressions illustrating the power-law behaviors 
in Eqs. |I| and ^ for the conductance G can be derived for the limit T <^ T^c = ft'^Bc/ks 
and V < Vrc = huRc/ e: 
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G{V,T)^^-^^ = GoiT)f(^^ 



(3) 



where Go{T) is the zero-bias conductance given by 



Rx \ hujc J r (a + 1) 
with ujci^ujRc) being the cutoff frequency of excitations, and 



(4) 



r 






r 


+ 1)' 


r 





(5) 



r(a;) is the Gamma function. Equation ^ exphcitly shows the power-law temperature de- 
pendence of Gq{= 1/Ro) for eV/2TTkBT <^ 1. For eV/2TTkBT ^ 1, the Gamma functions in 
Eq.|^ can be expanded, yielding the voltage power law GiV) oc V^"~^\ 

Equation H shows that, at low bias and fixed Risoi, the conductance can be scaled to a 
function of V/T as seen in Fig.^. The main figure shows this scaled conductance compared 
to the prediction of Eq.|^ (solid line). The arrows associated with the data sets at different 
temperatures indicate the points at which V = 0.25V^, where deviations from the low-voltage 
approximation begin to appear. For this scaling, Gq{T) has been used as a free parameter 
at each temperature. The resulting values for Gq{T) for both samples are compared with 
the numerical results (solid lines) in the inset of Fig. |] (Note that T is scaled to Trc for 
each sample.). Again, agreement is quite good. Sample 2, which has a larger Risoi (see 
Table j!]), shows the predicted decrease in Gq(T) and increase in the power-law exponent 
compared to sample 1. However, the low-temperature approximation (Eq. ^ dashed lines) 
deviates significantly from the data (and the finite-temperature theory) at T > 0.1 Trc. 
The scaling of voltage with temperature as well as the power-law dependence of / on are 
much more robust, with deviations of less than 15% from Eq.^ for V < 0.25V^, independent 
of temperature and Vjic- This fit, along with the strong temperature dependence of Go(T), 
indicates that the heating in Risoi is not significant, at least for V < 0.25V^. For T < 
60 mK, the data become nearly independent of T. We cannot presently rule out experimental 
artifacts as the explanation for this effect, which requires further study. 



In summary, we have measured the /-^characteristics of resistively isolated tunnel junc- 
tions with Rjsoi ^ Rk- We find strong Coulomb blockade with an increase of up to 10^ in 
the zero-bias resistance compared to the high-bias tunnel resistance of the junctions. For 
temperatures above 60 mK, the results agree well with the theory based on an "environ- 
mental" description of the resistor as a collection of oscillators with a linear spectrum. In 
particular we clearly observe the predicted power- law dependences of I{V) and Gq{T) with 
the power-law exponents dependent on Risoi in a manner consistent with the theory. 

We thank K.K. Likharev and A. Rylyakov for useful and stimulating discussions. This 
work was supported in part by AFOSR. 
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FIGURES 




v[Vd 

FIG. 1. I-V characteristics (thin sohd hne) of sample 1 at 70mK compared with the theoretical 
fit (thick dashed line) using Rjsol as a free parameter. The upper inset shows the schematic of 
the sample, where the shaded areas are the overlapping regions formed by Al shadow evaporation 
and the darkest areas indicate the junctions. C and Cs are the capacitances of the main junction 
and the secondary junction respectfully, and Risol represent the isolating resistors. The lower inset 
gives the linearly scaled / -V curve and associated fit. 
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V[nV] 

FIG. 2. I-V curves at different temperatures for sample 1(2) compared with numerical results 
from the full theory. The temperatures are 70, 90, 118, 135, 184 and 225mK (bottom to top) for 
sample 1 and 75, 90, 122, 143, 195and 225mK (bottom to top) for sample 2. 
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FIG. 3. Scaled conductance of sample 2 compared with the scaling function given by Eq. ^ 
(solid line) for data at 75, 90, 122, 143, 195and 225mK (top to bottom). The arrows indicate the 
points where V = 0.251c for the data set corresponding to the associated symbol. The inset shows 
the dependence of Go on temperature for both samples (squares for sample 1, circles for sample 2) 
compared with finite-temperature theory for each sample (solid lines). The dashed lines, which are 
shifted to asymptotically match the simulations, show the power-law behavior predicted by Eq. ^ 
for T^Trc. 
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TABLES 



TABLE I. Characteristics of the samples. 



O 1 AT 

bample No. 


1 


2 




170 


152 


F, (K\ 
^c\J^ ) 




4 8Q 




455 


422 


VRcitiV) 


71 


41 


Trc{K) 


0.84 


0.48 


Risoi{Trneas.){kQ.) 


40 


75 


Risol{fit){kn) 


52 


84 
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